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Comparison of dust-to-gas ratios in luminous, 
ultraluminous, and hyperluminous infrared galaxies 
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The dust-to-gas ratios in three different samples of luminous, ultraluminous, and hyperluminous infrared galaxies are 
calculated by modelling their radio to soft X-ray spectral energy distributions (SED) using composite models which 
account for the photoionizing radiation from HII regions, starbursts, or AGNs, and for shocks. The models are limited to a 
set which broadly reproduces the mid-IR fine structure line ratios of local, IR bright, starburst galaxies. The results show 
that two types of clouds contribute to the IR emission. Those characterized by low shock velocities and low preshock 
de nsities explain the far-IR dust emission, while those with higher velocities and densities contribute to the mid-IR dust 
emission. Clouds with shock velocities of 500 kms _1 prevail in hyperluminous infrared galaxies. An AGN is found in 
nearly all of the ultraluminous infrared galaxies and in half of the luminous infrared galaxies of the sample. High IR 
luminosities depend on dust-to-gas ratios as high as ~ 0.1 by mass, however most hyperluminous IR galaxies show dust- 
to-gas ratios much lower than those calculated for the luminous and ultraluminous IR galaxies. 
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1 Introduction 

Observations in the infrared (IR) have revealed galaxies with 
very high IR luminosities, up to Ljjj ~ 10 13 L Q (Lutz et 
al. 1998; Laurent et al. 2000; Tran et al. 2001; Verma et 
al. 2002). The origin of this high IR luminosity is under 
investigation (Genzel et al. 1998; Rowan-Robinson 1992, 
2000; Rigopoulou et al. 1999; Farrah et al. 2002; Verma et 
al. 2003; etc.). The main questions concern the type (such 
as HII galaxy, starburst nucleus galaxy, or active galactic 
nucleus (AGN)) of the IR bright and luminous galaxies, 
the source of dust illumination and heating by a photoion- 
izing flux and/or collisional processes, the composition of 
dust ranging from relatively large silicate grains to small 
graphite, PAH, ices, etc. (Spoon et al. 2002). 

Dust plays an important role in the understanding of the 
near and distant Universe (Schneider, Ferrara & Salvaterra 
2004). Formation of grains in the ISM is extremely ineffi- 
cient, the preferred sites being the atmospheres of evolved 
low mass (M < 8 M Q ) stars from where it is transported into 
the ISM through stellar winds (Whittet 1992). However, at 
high redshifts, e.g. z > 5, the evolutionary time-scales of 
low-mass stars (0.1 — 1 Gyr) start to be comparable with the 
age of the Universe. Type II supernovae are the only known 
dust sources with evolutionary time-scales shorter than the 
Hubble time. Dust enrichment of the ISM must have oc- 
curred primarily on considerably shorter time-scales in the 
ejecta of supernova explosions. 
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In the starburst environments and in the outflow regions 
of AGN, dust dominates the IR emission by reprocessing 
radiation from the primary source (stars and/or an active 
nucleus). Collisional processes which were found to be so 
important to explain the emission-line and continuum spec- 
tra (e.g. the strong high-ionization level line fluxes, the soft 
X-ray emission, strong heating of dust, radio synchrotron 
emission, etc.) from the nuclear and circumnuclear regions 
of AGNs, starbursts and HII regions (e.g. Contini & Vie- 
gas 2000) dominate the emitted radiation flux in mergers. 
Particularly, it was demonstrated that collisional heating of 
dust grains by gas in supersonic velocity regimes, leads to 
relatively high temperatures which could explain the IR flux 
between 3 and 1000 /im (the "IR bump") in the contin- 
uum spectral energy distribution (SED) of starburst galax- 
ies (Contini & Contini 2003) and AGNs (Contini, Viegas & 
Prieto 2004). 

Our aim is to test whether the IR luminosity is directly 
related to the amount of dust relative to gas, namely the 
dust-to-gas ratio (d/ g) which can be considered as a key pa- 
rameter in modelling galaxies revealing high dust formation 
regions and tracing stellar evolution. 

We investigate three luminosity classes commonly used 
for IR luminous galaxies: hyper-, ultra-, and luminous in- 
frared galaxies (HLIRGs, ULIRGs, and LIRGs, respectively) 
which correspond to > 10 13 , > 10 12 and > 10 11 h^ 2 L©, 
according to the definition of Sanders & Mirabel (1996). 
Three samples are chosen as representative of the different 
luminosity classes : Farrah et al. (2002), Rigopoulou et al 
(1999), and Verma et al. (2003). Only three galaxies from 
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the Verma et al. (2003) sample correspond to the definition 
of Sanders & Mirabel, the other showing lower IR lumi- 
nosities belong to the IR bright galaxies. However, in the 
following we will refer to them together with the luminous 
IR galaxies among the LIRGs. The samples are briefly pre- 
sented in Sect. 2. 

The galaxies in the different samples are chosen only 
because they are IR luminous and not by their type (e.g. 
starburst, HII regions, AGN, etc) which, viceversa, we will 
try to find out by modelling. Galaxies are generally classi- 
fied through the analysis of the SED of the continuum and 
of the absorption and emission bands of dust main com- 
ponents, e.g. silicates and PAHs (Spoon et al. 2002), on 
top of it. Diagnostic diagrams are generally used to distin- 
guish between AGN and star forming regions on the basis 
of PAH to continuum intensity ratios (Genzel et al. 1998, 
Peeters et al 2004), of the ratio of the continuum flux at dif- 
ferent wavelengths (Genzel & Cesarsky 2000), PAH to fine 
structure lines (FSL) (Genzel et al 1998) and FSL diagrams 
(Sturm et al 2002). However, recent observations of mid-IR 
AGN spectra by the Spitzer Space Telescope (Weedman et 
al. 2005) found no spectral parameters that unambiguously 
distinguish AGNs and starbursts based only on the slopes of 
the continuum spectra. 

The IR luminosity of galaxies depends on both dust re- 
processed radiation in the IR and bremsstrahlung from cool 
gas. Dust IR emission is usually analysed independently of 
bremsstrahlung which covers the whole radio - X-ray fre- 
quency range. Indeed, in the IR domain dust emission dom- 
inates, but the modelling of the SED requires the consistent 
calculation of gas and dust spectra (Sect. 3). Previous mod- 
elling of AGNs and starburst galaxies (Viegas, Contini & 
Contini 1999; Contini, Viegas & Prieto 2004) has shown 
that many different conditions coexist in each galaxy. The 
observed continuum and emission-line spectra account for 
all of them. Modelling the continuum SED by multi-cloud 
models, we will be able to determine the d/g ratios and the 
relative contribution of the different mechanisms which are 
at work in the IR luminous "active" galaxies : AGN, star- 
bursts, supernova remnants, etc. The results are presented 
in Sect. 4, discussed in Sect. 5 and summarized in Sect. 6. 

2 The samples 

The sample of Verma et al. (2003) contains starburst regions 
that are optically obscured, therefore similar to dense star 
forming regions in ULIRGs and HLIRGs. The IR luminosi- 
ties are between 10 9 and 10 12 L©, except for IC 342 corre- 
sponding to 3.74 10 7 L Q which will be used for comparison. 
We will check whether an active nucleus cannot be excluded 
in some of these starburst galaxies, leading to hybrid types. 

ULIRGs dominate the local luminosity function of galax- 
ies with L> 10 12 £© (Sanders & Mirabel 1996 and refer- 
ences therein). Most of their bolometric luminosity is emit- 
ted in the far-IR (8-1000 fim) (Rigopoulou et al. 1999) and 



Table 1 Characteristics of the galaxies 





type" 


d b 

(Mpc) 


logLm 
(L©) 


LIRGs : 








NGC 253 


HII 


3.27 


10.23 


IC 342 


HII 


0.45 


7.57 


II Zw 40 


HII 


10.5 


9.3 


M82 


HII 


2.7 


10.18 


NGC 3256 


HII 


37.95 


11.4 


NGC 3690 


HII 


41.64 


11.53 


NGC 4038/9 


HII 


21.9 


10.65 


NGC 4945 


HII,Sy2 


7.47 


10.66 


NGC 5236 


HII 


6.92 


10.11 


NGC 5253 


HII 


5.39 


9.21 


NGC 7552 


HII,Liner 


21.15 


10.83 


ULIRGs : 








Mrk 1014 


QSO 


652.3 


12.50 


UGC5101 


IrS 


160. 


11.94 


Mrk 231 


IrS 


168.8 


12.48 


Mrk 273 


IrS 


148. 


12.07 


Arp 220 


IrS 


72. 


12.11 


NGC 6240 


IrS 


96. 


11.78 


HLIRGs : 








F00235+1024 


nl 


2320. 


13.04 c 


07380-2342 


nl 


1160. 


13.34 c 


F10026+4949 


Syl 


4480. 


13.8L 


F 12509+3 122 


QSO 


3120. 


13.26 c 


13279+3401 


QSO 


1714. 


12.88 c 


14026+4341 


Syl 


1280. 


12.90 c 


F14218+3845 


QSO 


4840. 


13.06 c 


F16124+3241 


nl 


3381. 


13.02 c 


EJ 1640+41 


QSO 


4400. 


12.90 c 


18216+6418 


Syl 


1200. 


13.14 c 



a from Verma et al. (2003) for LIRGs, Rigopoulou et al. (1999) 
for ULIRGs, and from Farrah et al. (2002) for HLIRGs. 
b adopting H =75 kms^Mpc -1 

c 1-1000 fj,m luminosities obtained from the best-fitting 
combined starburst-AGN models (Farrah et al. 2002, table 3) 

they show more prominent FIR emission than either star- 
burst galaxies or AGNs (Peeters et al. 2004). Most ULIRGs 
are very rich in dust and gas, and they are interacting sys- 
tems with distorted morphologies. Their optical spectra mimic 
those of Seyfert galaxies (Rigopoulou et al. 1999, Lutz et al. 
1999). The MIR spectra of ULIRGs are more complicated 
than the spectra of other galaxies due to the presence of co- 
pious amount of dust in their nuclear regions (Genzel et al. 
1998), leading to strong absorption features which may dis- 
tort any emission component (Sturm et al. 2002). 

Infrared, millimeter and radio characteristics of ULIRGs 
are similar to those of starburst galaxies, while nuclear op- 
tical emission-line spectra are similar to those of Seyfert 
galaxies. A central AGN and circumnuclear star formation 
may both play a role in the FIR emission (Rigopoulou et 
al 1999). At intermediate luminosities, some ULIRG spec- 
tra are PAH dominated, while others show signs of a 8 fim 
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broad peak, indicating the presence of deeply dust-enshrouded 
ionizing sources (Rigopoulou et al. 1999). On the other hand, 
the spectra at high FIR luminosities bear close resemblance 
to AGN hot dust continua, showing little or no sign of PAH 
features (Clavel et al. 2000). Of the ultraluminous IRAS se- 
lected in Sect. 3, galaxies, 70-80% are predominantly pow- 
ered by recently formed massive stars, 20-30% by AGN 
(Genzel et al. 1998; Laurent et al. 2000). From multi-band 
HST imaging data Borne et al. (2001) find that ULIRGs 
show a significant population of very bright knots, which 
are most likely produced during the starburst/merger event. 

At least some HST observations (Farrah et al. 2002) 
show that a small fraction of HLIRGs are merging galaxies. 
Mergers between evolved galaxies trigger dust enshrouded 
starburst and AGN activity (Sanders & Mirabel 1996). Al- 
ternatively, HLIRGs may be very young or even "primeval" 
galaxies (Rowan-Robinson 2000, Farrah et al. 2002). Rowan- 
Robinson argues that the bulk of emission in HLIRGs is due 
to starburst activity, implying star formation rates > 1000 
M yr" 1 . Verma et al. (2002) and Farrah et al (2002) claim 
that both starburst and AGN activities are required to ex- 
plain the total IR emission, the average starburst fraction 
being 35%, within a range spanning from 80% starburst- 
dominated to 80% AGN-dominated. 

The types, the distances from Earth calculated adopting 
H = 75 km s _1 Mpc _1 , and the luminosities of the galaxies 
in the different samples are given in Table 1 . 

3 Modelling the spectra of starburst galaxies 

We start by modelling the IR emission-line spectra observed 
by Verma et al. (2003) for a sample of starburst galaxies 
with relatively high IR fluxes and select the models (Sect. 

3.1) which best explain the trend of the line ratios (Sect. 

3.2) , because line spectra are far more constraining than the 
SEDs. So we derive the grid of models suitable to model the 
SEDs (Sect. 3.3) by the method of Contini & Contini (2003) 
and Contini, Viegas & Prieto (2004) for different samples of 
luminous IR galaxies, and AGNs, respectively. Namely, gas 
and dust are heated and ionized consistently by photoioniza- 
tion and shocks (Sect. 3.4), so the ratio of dust reprocessed 
radiation to bremsstrahlung in the IR depends on dj g (Sect. 
3.5). 

3.1 The models 

We use composite models which account consistently for 
photoionization from an external source and for shocks, in a 
plane-parallel geometry (see Contini & Contini 2003, Con- 
tini et.al. 2004). Richtmyer-Meshkov instability between two 
fluids at different densities driven by shock waves and Kelvin- 
Helmholtz instability lead to perturbed circular material which 
is broken up into relatively small filaments. (Graham & Zhang 
2000, Contini & Formiggini 2001). The code SUMA (Vie- 
gas & Contini 1994, Contini & Contini 2003, Contini et. 



al. 2004, and references therein), is adopted for the calcula- 
tions of the spectra. SUMA includes shock hydrodynamics 
which implies a leading role of collisional processes, as well 
as excitation and destruction of dust grains. 

Other codes deal with dust emission in the IR, in partic- 
ular GRASIL (Silva et al. 1998) and MAPPINGS (Groves 
et al. 2006 and references therein). GRASIL predicts time 
dependent SEDs of galaxies from X-ray to radio, including 
the effects of age selective extinction of stellar populations. 
The d/g ratio, assumed to scale linearly with the metallic- 
ity of the residual gas, determines the optical depth of the 
clouds. MAPPINGS adopts stochastic quantum heating of 
dust and a power-law size distribution of grains with index 
arising from shattering. 

In our models the photoionizing flux from an AGN, Fh (in 
photons cm~ 2 s _1 eV _1 at 1 Ryd), corresponds to a power- 
law. For starburst galaxies, the radiation flux from a star 
cluster of age (t) is characterized by the ionization parame- 
ter U. For HII regions, the radiation flux is a black body cor- 
responding to the stellar colour temperature T». The shocks, 
characterized by the shock velocity, V s , the preshock den- 
sity, no, and the preshock magnetic field, Bo (= 10~ 4 gauss 
for all models), are created by collision of the ejected matter 
with the interstellar clouds. The abundance of He, C, N, O, 
Ne, Mg, Si, S, A, CI, and Fe relative to H and the dust-to-gas 
ratio, d/g are also input parameters. 

The clouds move outwards from the galaxy centre in 
the NLR of AGN (see e.g. Heckman et al. 1981), and in 
the neighbourhood of starburts, therefore, the photoioniz- 
ing flux and the shock act on the opposite edges of the 
clouds whose geometrical thickness, D, plays an important 
role in model results. The structure of the gas downstream 
is determined by the shock. Grains are sputtered depend- 
ing on V s and no (Contini et al. 2004). The sputtering rate 
determines the distribution of the grain radius throughout 
a cloud. An initial grain radius a gr = 0.2 fim is generally 
adopted. Small grains (< 0.2 /im) deriving from sputtering 
of larger grains contribute in the integrated continuum SED. 

It was found (Schneider et al. 2004) that for all but the 
smallest progenitor stars, silicate grains are the dominant 
dust compound. The characteristic grain sizes range between 
10~ 3 /im and a few x 0.1 /im, depending on the grain species 
and on the mass of the stellar progenitor. Very small grains 
are destroyed in regions of high UV radiation intensity (Rowan- 
Robinson 1992), so the maximum grain temperature which 
would account for the near-IR emission is limited by sub- 
limation. Most of the grain survive in the narrow line re- 
gion (NLR), where the the cloud velocities are ~ 300-1000 
kms -1 , while small grains (< 0.1 /im) are easily sputtered 
for shock velocities > 200 km s . 

PAHs are responsible for no more than 22% of the total 
3-1100 /im emission in normal galaxies (Dale et al. 2001). 
In the shock dominated starburst and AGN regions, small 
grains such as graphite and PAH are rapidly destroyed by 
sputtering, leading to an even lower percentage. Moreover, a 
strong radiation flux dilutes the PAHs which are hardly seen 
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Fig. 1 Top left: comparison of observed line ratios [SIV]/[SIII] vs. [NeIII]/[NeII] (asterisks) with models for HII regions + shocks 
(filled squares). Values of V s and n are indicated near each model. Top right: comparison of observed line ratios [ArIII]/[ArII] vs. 
[NeIII]/NeII] with models for HII regions + shocks. Same symbols as for top left diagram. Bottom left: comparison of observed line 
ratios [SIV]/[SIII] vs. [NeIII]/[NeII] (asterisks) with models for starbursts + shocks. Numbers refer to the starburst age in Myr. Open 
symbols represent geometrically thin clouds with D = 0.03 pc, filled symbols D = 3 pc; hexagons correspond to V s = 300 kms -1 , 
no = 300 cm -3 , triangles to V s = 200 kms -1 , no = 200 cm -3 . Bottom right: comparison of observed line ratios [SIV]/[SIII] vs. 
[NeIII]/[NeII] (asterisks) with models for AGN + shocks. Numbers represent log F h . filled triangles: V s =200 kms -1 , n = 200, 
D = 3 pc; empty triangles: V s =200 kms -1 , n = 200, D = 0.3 pc; black hexagons: V s =300 kms -1 , n = 300, D = 3 pc; black 
squares: V s =500 km s _1 , n = 300, D — 3 pc. All models except where indicated correspond to U=10. 



in AGN (Sturm et al. 2000, 2002 and Laurent et al. 2000). 
PAH, with radius <0.01 /im survive only in the extended 
"cirrus" regions of starburst galaxies (see Contini & Contini 
2003), and very seldom in the NLR of AGNs. Stochastic 
heating leads to temperature fluctuations around 20 K of 
a ffr =200 A grains, but small grains (a gr = 25 A) can reach 
50 K and even 1000 K for a sr =5 A(Draine 2003). However, 
their contribution to the IR is negligible because the cirrus 
region is located in the edges of the extended underlying 



galaxy, where d/g is even lower than in the ISM. So only 
silicate are accounted for in our models. 



3.2 The line ratios 

Verma et al. (2003) present ISO Short Wavelength Spec- 
trometer (SWS) spectra between 2.38 and 45 /im for 12 
starburst regions located in eleven galaxies with 10 9 L Q 
< Lir 10 12 L Q . Radio and infrared measurements demon- 
strate that some of the most active star-forming regions are 
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Table 2 Single-cloud models (top), the fluxes (in erg cm" 2 s -1 ) calculated in the different ranges (middle) and the maximum tem- 
perature of the grains (bottom); V s : the shock velocity, n : the preshock density, D : the geometrical thickness of the cloud, U : the 
ionization parameter, d/g : the dust-to-gas ratio, Fh : the ionization power-law (a=1.5) flux (in number of photons cm -2 s -1 eV -1 at 1 
Ryd), F radi0br : the radio flux calculated by bremsstrahlung, Fm br : the IR flux calculated by bremsstrahlung, F optbr : the optical flux 
calculated by bremsstrahlung, F UVbr ■ the UV flux calculated by bremsstrahlung, F IRd : the IR flux calculated by dust, T d (max) (K) : 
the maximum temperature of grains. T*=10 4 K is used for models ml-m5. 
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Fig. 2 The profile of the electron temperature, T c , of the electron density, N e , of the grain temperature, Td, and of the grain radius, 
a gr , throughout a cloud for models ml (left), m5 (middle), and m7 (right). 



optically obscured, with most of the bolometric luminosity 
emerging in the IR. 

We analyse the spectra observed by Verma et al. (2003, 
table 2) using the grids of models presented by Contini & 
Viegas (2001a,b). We focus on the ratios of MIR fine struc- 
ture line transitions relative to Ne, S, and Ar. These ele- 
ments, in fact, appear in at least two ionization levels. C, Si, 
and Fe etc., which are important in tracing dusty regions, as 
they can be trapped into dust grains and depleted from the 
gaseous phase, correspond to only one level. 

We have chosen line ratios belonging to the same ele- 
ment and corresponding to different ionization levels, which 
do not depend on the relative abundances, but only on the 
physical parameters. Indeed the relative abundances of the 
most coolant elements can influence the line ratios calcu- 
lated with a given model, changing the distribution of the 



ion fractional abundance downstream. However, we are deal- 
ing here with neon, sulfur, and argon which are not impor- 
tant coolants. Moreover, Verma et al. (2003) conclude that 
Ne and Ar are overabundant by no more than a factor of 
~3 relative to the solar ones, while S is underabundant by a 
similar factor (see also Contini & Contini 2003). We adopt 
cosmic abundances (Allen 1973). 

The [SIV]10.5/[SIII]18.7versus [NeIII]15.55/[NeII]12.8 
emission-line ratios observed in starburst galaxies are com- 
pared in Fig. 1 (top left) with models calculated using a 
black body radiation flux which corresponds to HII regions 
and shocks. The models which best reproduce the data were 
calculated (Contini & Viegas 2001b) with a stellar colour 
temperatures T*=10,000K, a ionization parameter U = 10, 
and a geometrical thickness of 3 pc. The relatively high U 
indicates that the emitting clouds are close to the stars. Only 
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Fig. 3 Comparison of the model SEDs : mO (long dash - dotted), 
ml (dotted lines), m2 (short dashed), m3 (dot - short dashed), m4 
(solid), m5 (long dashed), m6 (thin long dash - short dashed), and 
m7 (thick long dash - short dashed). 



Fig. 4 Comparison of continua calculated using models which 
account i) for a black-body flux with U = 10 (long-dashed lines 
for V s = 150 kms - , and solid lines for V s = 500 kms -1 ), and 
ii) for a power-law flux with log(Fh ) = 10 (short-dashed lines for 
V s = 150 km s- 1 and dotted lines for V s = 500 km s" 1 ) 



NGC 4945 is explained by lower fluxes (U ~ 1). The ob- 
served line ratios (asterisks) follow the trend of models with 
increasing n and V s which are indicated near each model. 
Higher V s in young starbursts, e.g. for II Zw 40, can be 
related to strong stellar winds produced by young massive 
stars (e.g. Wolf-Rayet stars). 

In Fig. 1 (bottom left), the [SIV] 10.5 /[SIII]18.7 ver- 
sus [Nelll] 15.55 /[Nell] 12.8 observed line ratios are com- 
pared with models (Contini & Viegas 2001b) which take 
into account black-body radiation from stars with a temper- 
ature distribution corresponding to that of a cluster with a 
certain age. The emitting gas is ionized and heated by radi- 
ation from the stars and by shocks. Only galaxies showing 
[SIV]/[SIII] > 0.15 are explained by starburst models. The 
emitting clouds are mosty geometrically thin (D < 0.03 
pc), and the best fitting models were calculated with U < 
0.1, indicating that the clouds are rather far from the radia- 
tion source, and/or that the sources are relatively weak. The 
starburst ages in Myr are indicated near each model. Galax- 
ies with line ratios > 0.1 correspond generally to young 
ages (t = 0.0 — 2.5 Myr), with a maximum of t = 3.3 
Myr for NGC 4038/4039. However, there is no clear trend 
between the observed line ratios and starburst ages or shock 
velocities. The shock velocities are generally of the order of 
200-300 kms- 1 . 

In Fig. 1 (bottom right) the galaxy sample is compared 
with models which account for radiation from an AGN and 
shocks. Also in these case the galaxies with low [SIV]/[SIII] 
cannot be explained by the models, while there is a rough 
trend of higher line ratios corresponding to higher flux in- 
tensity, Fh from the active centre. 



Fig. 1 (top left) shows that models calculated for HII 
regions explain the line ratios of all the sample galaxies. 
For sake of consistency, we have applied the same models 
to [ArIII]8.99 / [ArII]6.99 versus [Nelll] 15.5 /[Nell] 12.8 
diagrams (Fig. 1 top right). Indeed, the models fit the obser- 
vational data well. The distribution of the galaxies in the 
diagram is not very different from that corresponding to 
[SIV]/[SIII] versus [NeIII]/[NeII], and can be explained by 
the same shock velocity trend. 

The starburst galaxies observed by Verma et al. (2003) 
also appear in the sample of Lutz at al (1998). In a previous 
analysis by Viegas, Contini, & Contini (1999), the mod- 
elling of II Zw 40 and NGC 5253 led to gas with low V s 
and n , in agreement with their young age, on the basis of 
high [OIV]/([NeII]+0.44[NeIII]) versus [NeIII]/[NeII] line 
ratios. However, in the analysis of [SIV]/[SIII], we find that 
the same two galaxies are explained by HII region clouds 
with relatively high V s (500 km s _1 ), by clouds in starbursts 
with an age of ~ 2.5 Myr and velocities of 200-300 km s _1 , 
and by clouds in AGNs with velocities of ~ 300 km s _1 . 

Recall that ionization potentials of sulphur are lower 
than those of the corresponding levels of oxygen, therefore 
the stratification downstream of S +3 and S +2 ions is dif- 
ferent than that of +3 , Ne+ 2 , and Ne +1 . The ionization 
potentials of sulphur are also lower than those of the corre- 
sponding levels of argon, therefore the stratification down- 
stream of the S+ 3 and S +2 ions corresponds roughly to that 
of Ar +2 and Ar +1 . Therefore, the two top diagrams of Fig. 
1 show similar results, which are reasonably different from 
those obtained from the [OIV]/([NeII]+0.44[NeIII]) analy- 
sis. This indicates that clouds in different physical condi- 



© WILEY- VCH Verlag GmbH & Co. KGaA, Wcinhcim 



www.an-journal.org 



Astron. Nachr. / AN () 



7 



tions contribute to the different lines, e.g. [SIV] and [OIV], 
with various importance. Multi-cloud models were indeed 
suggested by Viegas et al. (1999, table 2) by modelling other 
galaxies of the sample, M 82, NGC 3256, and NGC 253, 
on the basis of infrared line ratios. Particularly, for M 82 
the contribution of the models corresponding to V s =400 
kms _1 and 100 kms _1 is 0.82 and 0.13 to the [OIVJ25.9 
line, but it is 0.40 and 0.45, respectively, to [Nelll] 15.5. 

Summarizing, the lines which appear in Fig. 1 diagrams 
are emitted by clouds located near the stars, which are char- 
acterized by temperatures of T* ~ 10 000 K. Model results 
show that H1I regions explain the [SIV]/[SIII] line ratios for 
all the galaxies of the sample, in agreement with the galaxy 
spectral types given by Verma et al (2003, Table 1). Clus- 
ter of stars located in different regions of the galaxies also 
contribute, as well as an active nucleus in some galaxies. 
More particularly, the galaxies of the Verma et al. (2003) 
sample show that II Zw 40, NGC 5253, NGC 4038/4039, 
NGC 3690A and NGC 3690 B/C present a multiple nature 
of HII regions, AGN, and starburst galaxies. 

3.3 The grid 

From the current literature and from diagnostic diagrams 
(Fig. 1) it became evident that many mechanisms coexist in 
each galaxy. Therefore, the grid of selected models (Table 
2) accounts for starbursts, AGN, HII regions, SNR winds, 
etc. We reduced the full set of models to a few prototypes to 
have clear results. 

It is generally believed that starbursts characterize lu- 
minous IR galaxies because dust is formed in the atmo- 
sphere of massive stars. So the models suitable to the star- 
burst clouds (m2-m5) show different physical parameters 
which were selected cross-checking the SEDs. An AGN, is 
often discovered in starburst galaxies, dominating the galac- 
tic central region where densities and velocities are rela- 
tively high, while the spectra emitted from shocked clouds 
in the external NLR of AGN is hardly distinguished from 
those emitted by clouds in the neighbourhood of starbursts. 
Therefore, we included in the grid only one model (m7) 
corresponding to the AGN, with large grains (a gr =1 jitm) 
which can survive sputtering downstream at relative high 
shock velocities and densities. 

SNR are strong emitters of radio and X-ray downstream 
of strong shocks created by the SN explosion. We include 
in the grid model m6, a shock dominated model (U=0) with 
a relatively high velocity suitable to explain soft-X-ray ob- 
servations. 

Finally, the underlying galaxy spectra must be taken into 
account, represented by models showing a lower V s , a low 
geometrical thickness, D due to the fragmentary charac- 
ter of the clouds, and a diluted photoionizing flux from the 
source. We have added to the grid model mO (d/g = 4 10~ 4 ) 
which was found necessary to reproduce the MIR contin- 
uum of most of the galaxies in the Alonso-Herrero et al. 
(2003) sample (Contini, Viegas & Prieto 2004). The d/g 
values in the underlying galaxy are typical of the ISM. 



The models (Contini & Viegas (2001a,b) adopted to ex- 
plain the line ratios, were calculated assuming dust-to-gas 
ratios of 4 10 -5 by mass, which are low, even compared to 
the Galactic ISM (d/g ~4x 10~ 4 ). By combining these 
models to fit the continuum SED, more realistic d/g ratios 
will be obtained. 

The comparison of calculated with observed SEDs indi- 
cates that the best agreement is obtained by d/ g higher by 
some orders of magnitude in all the models, except in those 
corresponding to low V s (100-150 kms -1 ) and low no (100 
cm -3 ). The models with V s > 200 km ^recalculated adopt- 
ing a more suitable d/g however, could not explain the ob- 
served line ratios (Figs. 1), indicating that the line spectra 
are emitted from gas less obscured by dust and that dust is 
not distributed homogeneously throughout a galaxy. Lines 
emitted from the dusty regions are weak compared to those 
emitted from dust-free ones because a high d/g speeding 
up the cooling rate, reduces the ion fractional abundances 
by recombination. 

3.4 The single-cloud models 

The calculated bremsstrahlung integrated throughout the ra- 
dio, IR, optical, U V, and X-ray ranges for the different mod- 
els, are given in rows 7, 9, 10, 11, and 12 of Table 2, re- 
spectively. In row 9 the calculated flux emitted in the IR 
by reprocessed dust, FiR d is shown for comparison. Notice 
that both bremsstrahlung and reradiation by dust increase 
with V s , no and U. Particularly, shock-dominated mod- 
els (calculated with U = 0) such as m6 show very low 
bremsstrahlung radiation. 

The maximum temperatures of the grains downstream 
are given in the last row of Table 2 (bottom). The tempera- 
tures increase with V s . However, density and grain sputter- 
ing also affect the heating and cooling of dust. This can be 
noticed by comparing the maximum temperature of grains 
calculated with models mO and ml. 

The profiles of the temperature and density of the gas, 
and of the temperature and radius of dust grains calculated 
for models ml, m5, and m7 are shown in Fig. 2 (left, middle, 
and right diagrams, respectively) for comparison. 

The left edge of the diagrams represents the shock front. 
Here collisional processes dominate the heating of dust and 
gas up to a certain distance downstream which is determined 
by the shock velocity and by the density. The grains are de- 
stroyed by sputtering in the region close to the shock front 
(middle and right diagrams). The photoioning flux reaches 
the opposite edge, heating and ionizing gas and dust. 

The continuum SEDs calculated by models m0-m7 are 
compared in Fig. 3. Two curves correspond to each model, 
one representing dust reradiation and the other bremsstrahlung. 
The bremsstrahlung peak between 10 13 and 10 15 Hz de- 
pends on photoionization, while the maximum at higher fre- 
quencies is determined by the shock. For each model, the 
peak intensity of dust reprocessed radiation depends on d/ g 
and the peak frequency on V s . However, Fig. 3 shows that 
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Fig. 5 Modelling the SEDs of the Verma et al. (2003) luminous IR galaxies sample. Black triangles: data from the NED. Black squares 
in the II Zw 40 diagram : data from Galliano (2004). Dotted lines: ml; short-dashed lines: m2-m3; solid lines: m4; long-dashed lines: 
m5; long dash-dot lines: m7; short-long-dashed lines: m6; dash-dotted line: bb model for the old star population; solid line: synchrotron 
radiation. 



© WILEY- VCH Verlag GmbH & Co. KGaA, Wcinhcim 



www.an-journal.org 



Astron. Nachr. / AN () 



9 




log(^) log(f) log(f) 

Fig. 6 Modelling the SED of the ULIRGs from the Rigopoulou et al. (1999) sample. Symbols as in Fig. 5. Empty triangles: Rigopoulou 
et al. (1999) IR data at 5.9 and 7.7 ^m 



sputtering may change the distribution of the grains and literature to fulfil the dataset of ULIRGs and HLIRGs by 
shift the peak at lower frequencies (e.g. model m7). The Rigopoulou et al and Farrah et al., respectively. The data 
MIR flux emitted by dust corresponding to V s = 300 km s -1 andcorrespond to different apertures. In the diagrams we will 



500 kms -1 (models m4 and m5, respectively) show very 
similar SEDs, while the corresponding bremsstrahlung at 
frequencies > 10 13 Hz are different. Therefore, we will con- 
strain the models by both the IR dust reradiation flux and the 
bremsstrahlung in the IR, optical, UV, and X-ray range. 

The bremsstrahlung of models m6 and m7, both cal- 
culated with V s =1000 kms -1 , show different distributions 
because m6 is shock dominated while m7 is radiation dom- 
inated. Moreover, they are calculated by different n and 
different grain sizes. 

The accuracy of the data in the optical-UV range is not 
high enough to select between power-law and black-body 
fluxes by modelling the continuum SED only (see also Con- 
tini & Viegas 2000, Fig. 1). Shock velocities instead can be 
estimated by the peak at high frequencies of bremsstrahlung 
and by the frequency of the IR emission peak (Contini et al. 
2004). A comparison of different models is shown in Fig. 4. 

3.5 The continuum SED 
3.5.1 Selection of the data 

The data from the NED (NASA Extragalactic Database) 
are adopted for LIRGs, complemented with data from the 



distinguish between those taken with the smallest aperture 
(up to 14 arcseconds, black triangles), those with large aper- 
tures (white squares), and those integrated from maps (black 
squares). The observational errors are not shown in the dia- 
grams in order to obtain less confused pictures. 

All the data are important in modelling. Those taken 
with a large aperture cover different regions characterised 
by radiation from different sources in the same galaxy, there- 
fore they correspond to multi-cloud models. This method 
led for example to the identification of the starbursts in the 
region close to the nucleus of the low luminosity AGN NGC 
4579 (Contini 2004), or in the outskirts of NGC 7130 AGN 
(Contini et al. 2002). Data obtained with different apertures 
at frequencies close to silicate, PAH, graphite , ice, etc. bands 
may include different absorption/emission fluxes (Angeloni 
et al. 2007, in preparation). 

We try to reproduce ("by eye") the largest number of 
data, constraining the modelling by upper limits. The curves 
representing dust emission and bremsstrahlung from differ- 
ent clouds are not summed up in the diagrams in order to 
compare directly each model with the dataset trend. 

The old stellar population background is represented by 
a black body (bb) with T~ 3000-5000 K. The data contam- 
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Fig. 7 Modelling the SED of the HLIRGs from the Farrah et al. (2002) sample (all fluxes at rest-frame wavelengths shortward of 4 
/Ltm are treated as upper limits due to unqualified contribution from the host galaxy). Symbols as in Fig. 5. The filled circle represents 
the upper limit datum in the X-ray range by Wilman et al. (1998). 
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inated by different contributions of light from the stars, are 
nested within the bb curve, e.g. NGC 21 10 (Contini & Vie- 
gas 2000), Arp 220 (Contini et al. 2004). The less contam- 
inated data show the very bremsstrahlung from the clouds 
and constrain the models. 

For all the galaxies, X-rays are upper limits because 
the fluxes are integrated on a large frequency range. The 
power-law trend of the data in the radio continuum of many 
galaxies is explained by synchrotron radiation created by 
the Fermi mechanism at the shock front. 

The IR dust reradiation "bump" is generally complex, 
accounting for the contribution of different grains in differ- 
ent conditions. 

Finally, the IR bands on top of the continuum are not 
included in the calculations of dust, therefore, they can be 
recognised comparing calculated with observed SEDs. 

3.5.2 Modelling approach 

Models are calculated at the nebula, while data are observed 
at Earth. Therefore to suite model results to the observed 
SED, the results must be multiplied by a factor jj which ac- 
counts for the distance to Earth of the galaxy, for the dis- 
tance of clouds from the galaxy centre, and for the filling 
factor. 

To determine the dust-to-gas ratios in each cloud of a 
galaxy, we will use the method presented for the starburst 
regions in luminous infrared galaxies by Contini & Contini 
(2003). Dust and gas are coupled throughout the shock front 
and downstream in each cloud and the IR flux from dust is 
consistently calculated with bremsstrahlung, so the 77 factors 
corresponding to IR emission by dust and to bremsstrahlung 



are included in the galaxy field and corrections have not 
been applied (see e.g. Code & Welch, 1982). 

A deep feature at ~ 10 fira indicates strong absorption 
by silicates, particularly for in NGC 253, IC 342, M 82, 
NGC 5236, and NGC 5253. For NGC 5236 and NGC 7552 
there is also evidence of PAH emission in the mid-IR. 

The wide peak in the IR, observed in all the galaxies 
except NGC 4038/4039, can be explained by the sum of at 
least two models corresponding to different shock veloci- 
ties. The emission which peaks longward of 95 fim cor- 
responds to low-velocity shocks (V s = 100-200 kms" 1 ), 
which leads to maximum grain temperatures of ~ 50 K, 
while the emission which peaks at ~ 30 /jm corresponds to 
relatively high velocity shocks (> 300 kms -1 ) with maxi- 
mum grain temperatures of ~ 120 K. We have tried to obtain 
an agreement with the data within the observational error 
(< 30 %). For some galaxies the data are insufficient to con- 
strain the models and clouds with V s = 500 km s~ 1 cannot be 
excluded. However, they correspond to lower relative con- 
tributions. These results are consistent with those found for 
starbursts and LIRGs by Contini & Contini (2003). 

Fig. 5 shows that the IR reradiation by dust of NGC 253, 
IC 342, NGC 3256, NGC 3690, NGC 5236, and NGC 7552 
can be also modelled by a unique model corresponding to 
a AGN with V s = 1000 kms" 1 , n = 1000 cni~ 3 , and a 
power-law flux with logFh =11 (model m7) as was found 
in Arp 220 (Contini et al. 2004, Fig. Al). Indeed, the IR 
emission extends all over the IR wavelength range because 
the grains are heated to high temperatures in the immediate 
post-shock region. The grains are large enough (a gr = 1 /jm) 
to survive sputtering, and they are heated to temperatures of 
~ 50 K by the radiation flux from the active center on the 



must be equal. The is determined by fitting the bremsstrahlung^PP° site side of the cloud ^ 2 >' Consider for instance 

u- l * j *■ j »v t t u T4GC 253, besides the good agreement to the small aperture 

which extends from radio to X-rays. To recover the same 77 



for the IR bump, we change the d/g input parameter, which 
affects directly dust reradiation. These changes, however, 
must account for the multi-cloud model. 

This procedure is straightforward, but not trivial because 
changing the d/g ratio by a factor > 10 leads to a different 
SED of the bremsstrahlung. Therefore the d/g ratios result 
from cross-checking the dust reprocessed radiation by the 
bremsstrahlung, for each cloud in every galaxy. 

4 Results 

The results of modelling of the continuum SEDs of LIRGs, 
ULIRGs, and of HLIRG are presented in the following sec- 
tions. For most objects one type of clouds is not enough to 
explain the whole dataset from FIR to near-IR. 

4.1 LIRGs 

The modelling of LIRG is shown in Fig. 5. The data (see 
Appendix) come from the NED. These data are generally 
corrected for galactic reddening, however, background stars 



aperture 

data by the AGN model in the optical range and by the low 
velocity model (ml) in the UV, we must account also for the 
data observed with large apertures which are explained by a 
starburst model (m4). 

Synchrotron emission in the radio is observed in all the 
galaxies, except IC 342, II ZW 40, and NGC 5253. Free- 
free continuum absorption is evident in a few objects (NGC 
3256, NGC 3690, NGC 7552). 

The hybrid nature of LIRGs (starburst+AGN) appears in 
about half of the objects, (see Table 3). 

4.2 ULIRGs 

Rigopoulouetal. (1999) present, in their Table 1JSOPHOT- 
S spectroscopy for a sample of ULIRGs. We have selected 
those which have a rich number of data from the NED (see 
Appendix), and we have added the two data points at 5.9 
and 7.7 fim reported by Rigopoulou et al. (1999). 

In the absence of MIR spectroscopy for ULIRGs and 
HLIRGs we assume the same grid of models used for LIRGs 
are suitable for modelling systems of higher luminosity. The 
results are shown in Fig. 6. 
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Table 3 Log(?7) a and the relative contribution of the different clouds, W; for LIRGs 
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0.11 



u.u 



u.u 


8^ 


TTV 
u V 


n n 


n in 
U. 1U 



u.u 


ft ft 

u.u 


U.UJU 



u.u 



u.u 


8S 

U.O J 


X-ray 


0.0 


0.477 


0.0 


0.0 


0.0159 


0.0 


0.460 


0.0477 


IR d 


0.0 


0.21 


0.0 


0.0 


0.41 


0.0 


0.0 


0.37 


NGC5253 




-5.1 






-11.7 




-13.2 




opt 


0.0 


0.607 


0.0 


0.0 


0.39 


0.0 


0.0 


0.0 


UV 


0.0 


0.87 


0.0 


0.0 


0.12 


0.0 


0.0 


0.0 


X-ray 


0.0 


0.860 


0.0 


0.0 


0.00904 


0.0 


0.131 


0.0 


IRd 


0.0 


0.62 


0.0 


0.0 


0.37 


0.0 


0.0 


0.0 


NGC7552 








-11.5 


-10.8 




-13. 


-12.7 


opt 


0.0 


0.0 


0.0 


0.011 


0.29 


0.0 


0.0 


0.69 


UV 


0.0 


0.0 


0.0 


0.0 


0.15 


0.0 


0.0 


0.84 


X-ray 


0.0 


0.0 


0.0 


0.0059 


0.205 


0.0 


0.594 


0.195 


IRd 


0.0 


0.0 


0.0 


0.47 


0.41 


0.0 


0.0 


0.11 



a in the first row of each galaxy 
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Table 4 Log(?7) a and the relative contribution of the different clouds, W; for ULIRGs 





under 


-gal 


starbursts 


winds 


AGN 




mO 


ml 


ml 


m3 


m4 


m5 


m6 


m7 


MRK1014 


- 


-7. 


_ 


_ 


_ 


-12.9 


-12. 


- 


opt 


0.0 


0.16 


0.0 


0.0 


0.0 


0.83 


0.0044 


0.0 


UV 


0.0 


0.33 


0.0 


0.0 


0.0 


0.65 


0.012 


0.0 


X-ray 


0.0 


0.0051 


0.0 


0.0 


0.0 


0.0149 


0.98 


0.0 


IR d 


0.0 


0.019 


0.0 


0.0 


0.0 


0.98 


0.0 


0.0 


UCG5101 


-9.1 


- 


. 


. 


. 


-13.3 


- 


-13.7 


opt 


0.0033 


0.0 


0.0 


0.0 


0.0 


0.02 


0.0 


0.97 


UV 


0.0047 


0.0 


0.0 


0.0 


0.0 


0.015 


0.0 


0.98 


X-ray 


0.0 


0.0 


0.0 


0.0 


0.0 


0.649 


0.0 


0.351 


IRd 


0.0078 


0.0 


0.0 


0.0 


0.0 


0.64 


0.0 


0.35 


MRK231 


-8. 


- 


_ 


_ 


_ 


-11.8 


- 


-13.4 


opt 


0.016 


0.0 


0.0 


0.0 


0.0 


0.24 


0.0 


0.74 


UV 


0.024 


0.0 


0.0 


0.0 


0.0 


0.19 


0.0 


0.78 


X-ray 


0.0 


0.0 


0.0 


0.0 


0.0 


0.967 


0.0 


0.033 


IRd 


0.0046 


0.0 


0.0 


0.0 


0.0 


0.96 


0.0 


0.033 


MRK273 


- 


- 


_ 


_ 


_ 


-12.9 


- 


-13.5 


opt 


0.0 


0.0 


0.0 


0.0 


0.0 


0.031 


0.0 


0.97 


UV 


0.0 


0.0 


0.0 


0.0 


0.0 


0.024 


0.0 


0.97 


X-ray 


0.0 


0.0 


0.0 


0.0 


0.0 


0.745 


0.0 


0.255 


IRd 


0.0 


0.0 


0.0 


0.0 


0.0 


0.74 


0.0 


0.25 


Arp 220 


- 


-9.1 








-13.3 


- 


-12.9 


opt 


0.0 


0.0001 


0.0 


0.0 


0.0 


0.007 


0.0 


0.99 


UV 


0.0 


0.0003 


0.0 


0.0 


0.0 


0.003 


0.0 


0.997 


X-ray 


0.0 


0.0 


0.0 


0.0 


0.0 


0.23 


0.0 


0.77 


IRd 


0.0 


0.001 


0.0 


0.0 


0.0 


0.22 


0.0 


0.78 


NCG6240 














-12.4 


-13. 


opt 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


UV 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


1.0 


X-ray 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.960 


0.0397 


IRd 


0.0 


0.0 


0.0 


0.0 


0.0 


0.0 


0.15 


0.85 



a in the first row of each galaxy 



For all the galaxies of the ULIRG sample, except Mrk 
1014, the low d/g clouds (ml) which significantly contribute 
in the LIRGs are not found, while an AGN nucleus is re- 
vealed by modelling the IR prominent flux. Free-free con- 
tinuum self-absorption is evident in NGC 1014 and Mrk 
231. 



Robinson (1992). Farrah et al. (2002) found that AGN and 
starburst luminosities correlate, suggesting that a common 
physical factor, plausibly the dust mass, governs the lumi- 
nosities of starbursts and AGNs in HLIRGs. 



4.3 HLIRGs 

At the brightest end of the luminous IR galaxy population 
lie the HLIRGs, with Ltr > 10 13 h^ 2 L Q (Rowan-Robinson 
2000; Farrah et al. 2002; Sanders & Mirabel 1996; Tran et 
al. 2001; Verma et al. 2002). They claim that these objects 
may reveal an entirely new class of IR-emitter objects. 

Farrah et al. (2002) present sub-millimetre photometry 
for eleven HLIRGs and use models by Efstathiou, Rowan- 
Robinson, & Siebenmorgen (2000) for starbursts and by Ef- 
stathiou & Rowan-Robinson (1995) for AGNs to examine 
the nature of the IR emission. Dust composition for star- 
bursts is given by the dust grain model of Siebenmorgen & 
Krugel (1992) and by the multigrain dust model of Rowan- 



The models we use here are different due to the impor- 
tant role of collisional processes. The results of our mod- 
elling are shown in Fig. 7. We found that the IR flux can 
be explained by two types of clouds, one corresponding to 
low V s - n values, low U, and low d/g, (ml) while the 
second type of clouds corresponds to rather high velocities 
V s = 500 kms^ 1 , high d/g values, and a high ionization 
parameter U, (m5). The models correspond to Td(max) of 
~ 50 K and 120 K, respectively. The X-ray upper limit for 
the HLIRG F00235+1024 (Wiman et al 1998) integrated in 
the 0.1-2.4 keV band, is well explained by model m5 calcu- 
lated by V s =500 km s _1 . Free-free continuum absorption is 
evident in P070380-2342 and P18216+6418. 
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Table 5 Log(?7) a and the IR d relative contribution of the different clouds, Wj for HLIRGs 





under-gal 


starbursts 


winds 


AGN 




mO 


ml 


ml 


m3 


m4 


m5 


m6 


ml 


F00235+1024 


_ 


-7.9 


_ 


_ 


_ 


-14.6 


- 


- 


IR d 


0. 


0.11 


0. 


0. 


0. 


0.89 


0. 


0. 


07380-2342 


_ 


-7.9 


_ 


_ 


_ 


-13.7 


- 


- 


IRd 


0. 


0.01 


0. 


0. 


0. 


0.98 


0. 


0. 


F10026+4949 


_ 


-8. 


_ 


_ 


_ 


-14.7 


- 


- 


IRd 


0. 


0.11 


0. 


0. 


0. 


0.89 


0. 


0. 


F12509+3122 


_ 


-6.7 


_ 


_ 


_ 


-15. 


- 


- 


IRd 


0. 


0.84 


0. 


0. 


0. 


0.16 


0. 


0. 


13279+3401 


_ 


-6.8 


_ 


_ 


_ 


-12.8 


- 


- 


IRd 


0. 


0.05 


0. 


0. 


0. 


0.95 


0. 


0. 


14026+4341 




-6.5 








-13.5 


- 


- 


IRd 


0. 


0.21 


0. 


0. 


0. 


0.79 


0. 


0. 


F14218+3845 




-8. 








-14.6 






IRd 


0. 


0.09 


0. 


0. 


0. 


0.91 


0. 


0. 


F16124+3241 




-8.4 








-14.6 






IRd 


0. 


0.04 


0. 


0. 


0. 


0.96 


0. 


0. 


EJ1640+41 




-8.4 








-14.4 






IRd 


0. 


0.02 


0. 


0. 


0. 


0.98 


0. 


0. 


18216+6418 




-6. 








-12.5 






IRd 


0. 


0.07 


0. 


0. 


0. 


0.93 


0. 


0. 



a in the first row of each galaxy 



5 Discussion 

5.1 The relative importance of the different energy 
sources in each galaxy 

The r\ factors (Sect. 3.5.2) derived by modelling the SED, 
allow to calculate the relative contribution of AGNs, star- 
bursts, high velocity winds, and the underlying galaxy, in 
the different frequency ranges A„ (radiobr, IRt, r , optbr, 
UVbr, Xbr, IRd, Table 2, bottom) for each galaxy, by W, = 
(Fa v )i £ (Fa v )i (v)i ■ They are given in Tables 3 and 4 
forLIRGs and ULIRGs, respectively. The IR bremsstrahlung 
is lower by a large factor than reradiation from dust, there- 
fore is neglected. 

The IRd contribution from the starburst clouds with dif- 
ferent velocity is generally high, except for M82 which shows 
a large relative contribution from the background galaxy. 
The strong IR,j from the starburst dust relative to AGN de- 
pends on the different character of the flux from the external 
source: a power-law in the AGN and a multiple black body 
in the starburts. Tables 3 and 4 show, on the other hand, 
that when an AGN is present, its contribution prevails in the 
optical and UV range, while in the X-ray the high velocity 
winds dominate. 

Tables 3 shows that relatively low V s - n clouds (V s 
< 200 kms" 1 ) dominate in M82 and NGC 5253. Starburst 
clouds with V s > 300 kins" 1 are less indicated in LIRGs. 
An AGN contributes to the IR by ~ 23 % in NGC 253 and 
IC 342, by ~ 40% in NGC 3256, NGC 3690, and NGC 
5236. In ULIRGs (Table 4) the AGN is always present (ex- 



cept for Mrk 1014) and dominates by > 75 % in Arp 220 
and NGC 6240. 

In Table 5 the relative contribution to the IR fluxes are 
given for HLIRGs. Most of the data in Fig. 6 are upper lim- 
its, so we focus on the IR. HLIRGs are modelled by two 
types of clouds. The relative contribution to IR^ of clouds 
with V s =500 kms" 1 (m5) is > 90% for all the galaxies of 
the sample, except for F12509+3122 which shows a large 
contribution from the underlying galaxy (ml). 

5.2 The distance of the emitting clouds from the 
galaxy center 

We have found by modelling the SEDs that the 77 factors cor- 
responding to models mO and ml are higher than those cor- 
responding to models m2-m7 by factors of 10 6 — 10 7 . This 
can be explained by larger distances of the clouds from the 
galaxy centre, and/or by a large number (N <~ l/ff where 
ff is the filling factor) of low V s - no clouds. The expla- 
nations are both valid. Indeed, low V s -no clouds (models 
mO and ml) with d/g ratios even lower than those of the 
Galactic ISM (Table 2) are characteristic of the extended 
underlying galaxy and low V s (< 200 kms -1 ) - n (~ 100 
cm~ 3 ) clouds were found in the outer narrow line region of 
AGNs which may also contain starburst regions (Contini et 
al. 2002). 

We define ?7=(d c /d) 2 (1/ 'ff), where d is the distance of 
the galaxy to Earth and d c the distance of clouds from the 
galaxy center. Fragmentation by shocks leads to clouds with 
different geometrical thickness, D (Table 2), particularly 
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Table 6 The distances of the clouds (in pc) from the center of each galaxy 





! under-gal 




starbursts 




| winds 


| AGN 


| radio | 




mO ml 


m2 


m3 m4 


m5 


m6 


ml 


m s ^ n 


LIRGs : 
















NGC 253 


220 


- 


20.4 




3.63 


2.9 


141.3 


IC 342 


- 


1.4 


1.26 


- 


0.56 


0.28 


6.3 


II Zw 40 


52 


- 


6.6 


- 


- 


- 


- 


M82 


85 


- 


10.7 


- 


4.8 


- 


141. 


NGC 3256 


370 


- 


85.1 


- 


30.2 


21.4 


631. 


NGC 3690 


420 


- 


100. 


- 


13.2 


16.6 


933. 


NGC 4038/4039 


- 


- 


- 


- 


12.3 


- 


275.4 


NGC 4945 


- 


59.6 


33.1 


- 


5.95 


- 


331. 


NGC 5236 


270 


- 


24.5 


- 


6.16 


4.9 


173.7 


NGC 5253 


150 


- 


7.58 


- 


1.35 


- 


- 


NGC 7552 


- 


- 


37.6 84.1 


- 


6.68 


9.44 


237.1 


ULIRGs : 
















Mrk 1014 


2000. 


- 


- 


229. 


653. 


58.2 


518.8 


UGC5101 


45 


- 


- 


35.5 


- 


22.38 


- 


Mrk 231 


170 - 


- 


- 


211. 


- 


33.5 


1.7e3 


Mrk 273 


- 


- 


- 


52.5 


- 


26.3 


1.05e3 


Arp 220 


2029. 


- 


- 


16.11 


- 


25. 


- 


NGC 6240 


- 


- 


- 


- 


60.25 


30.2 


- 


HLIRGs : 
















F00235+1024 


2600 


- 


- 


113.5 


- 


- 


- 


07380-2342 


1300 






164. 








F10026+4949 


4500 






199. 








F12509+3122 


14000 






100. 








13279+3401 


6760 






676. 








14026+4341 


7200 






229. 








F1421 8+3845 


4800 






245. 








F16124+3241 


2100 






173.8 








EJ 1640+41 


2780 






278. 








18216+6418 


11800 






668. 









low in the underlying galaxy. We adopt ff~ D m i n ID max 
= 10~ 4 (D m i n and D max are the minimum and maximum 
geometrical thickness of the clouds in Table 2) for the un- 
derlying galaxy clouds (mO and ml), and j(f ~1 for the star- 
burst and AGN clouds (m2-m7). 

Models m6 and m7 with a high shock velocity (V s =1000 
kms -1 ) explain the observed soft X-rays, which are upper 
limits. Dust emission calculated by model m7, which ap- 
pears in most of the ULIRGs, explains the mid- to far-IR 
flux observed in some galaxies. The rj of models m6 and 
m7 are relatively low indicating that the high V s clouds are 
located closer to the centre. 

The distances of the clouds from the galaxy center are 
shown in Table 6. We have found by fitting the continuum 
SED of LIRGs that a large number of low V s , no clouds 
are located in the outskirt of the galaxies. Even adopting the 
low ff indicated above, the distance of these clouds from 
the center results higher by a factor > 10 than that of clouds 
corresponding to higher velocities and densities. The large 
regions of low V s - n clouds with low d/g are less evi- 
dent in ULIRGs, as high and low velocity clouds coexist be- 
cause merging leads to fragmentation and mixing. Indeed, 



the galaxies selected from the Rigopoulou et al. (1999) sam- 
ple are few but they all show an irregular morphology due 
to merging (e.g. Leech et al 1994). Table 6 shows that the 
average distances of the starburst clouds from the centre is 
similar for LIRGs and ULIRGs, within 100 pc for almost all 
galaxies, except for Mrk 1014. The luminous IR galaxies of 
the Verma et al sample, with L IR > 10 11 L©, NGC 3256 
and NGC 3690 which are identified as a colliding pair and 
a peculiar merger, respectively, show starburst cloud dis- 
tances from the centre ~ 100 pc. There are not enough data 
to model the NCG 4038/4039 interacting complex. 

HLIRG redshifts (Farrah et al. 2002, table 1 and ref- 
erences therein) are higher than those of LIRGs, leading to 
distances larger by factors of 100- 1000. We have found (Ta- 
ble 6, bottom) that d c should be higher for HLIRGs than for 
LIRGs by about a factor > 10 in average, adopting the same 
ff which were used for the other galaxy samples. To recover 
the dimensions of the extended region of the other galaxies, 
lower ff should be adopted. It seems that Mrk 1014 with log 
(L/^/Lq) = 12.5 has characteristics similar to HLIRGs. 

Regions of 200-900 pc are generally characteristic of 
the NLR of objects such as QSOs, Seyfert 1, and narrow 
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Table 7 Dust-to-gas ratios (in units of 4 x 10 ) in each galaxy 





under-gal 




starbursts 




| winds 


| AGN | 




mO ml 


m2 


m3 m4 


m5 


m6 


ml 


LIRGs : 














NGC 253 


0.70 


- 


100 


- 


100 


200 


IC 342 


- 


100 


100 


- 


100 


250 


II Zw 40 


0.75 


- 


158 


- 


- 


- 


M82 


6. 


- 


340 


- 


100 


- 


NGC 3256 


4.7 


- 


200 


- 


100 


63. 


NGC 3690 


3 


- 


100 


- 


100 


79. 


NGC 4038/9 


- 


- 


- 


- 


100 


- 


NGC 4945 


- 


100 


100 


- 


100 


- 


NGC 5236 


0.75 


- 


100 


- 


100 


126. 


NGC 5253 


0.24 


- 


100 


- 


100 


- 


NGC 7552 


- 


- 


8. 20. 


- 


100 


79. 


ULIRGs : 














Mrk 1014 


0.95 


- 


- 


5. 


- 


- 


UGC5101 


6. 


- 


- 


30. 


- 


100. 


Mrk 231 


0.16 - 


- 


- 


4. 


- 


126. 


Mrk 273 


- 


- 


- 


5. 


- 


100 


Arp 220 


0.15 


- 


- 


50. 


- 


100 


NGC 6240 


- 


- 


- 


- 


20 


31.6 


HLIRGs : 














F00235+1024 


2.4 


- 


- 


5. 


- 


- 


07380-2342 


7.5 






50. 






F 10026+4949 


1.2 






100. 






F12509+3122 


0.06 






100. 






13279+3401 


0.19 






0.6 






14026+4341 


0.06 






6.3 






F14218+3845 


0.95 






4. 






F16124+3241 


1.9 






16. 






EJ 1640+41 


0.95 






5. 






18216+6418 


0.04 






0.8 







a adopting H =75 km s *Mpc 1 
1-1000 fim luminosities obtained from the best-fitting combined starburst-AGN models (Farrah et al. 2002, table 3) 



line galaxies. This indicates that also HII regions are located 
within the NLR of these galaxies. Synchrotron emission in 
the radio range corresponds to larger distances from the cen- 
tre, which is reasonable enough considering that it is gener- 
ally created in the lobes of extended jets. 

5.3 The dust-to-gas ratios 

The d/g ratios relative to those of the Galactic ISM (<~ 
4 x 10~ 4 by mass) for each type of clouds in each galaxy are 
given in Table 7. More particularly, d/g as high as 0. 1 were 
found in M 82, and d/g > 0.04 in the starburst clouds of 
most of the LIRGs. This may indicate that dust in the star- 
burst regions of LIRGs is of SN origin. Indeed high velocity 
winds were found modelling the X-ray in Figs. 5. SN cor- 
responding to explosions of 10-30 M Q stars, are potentially 
the most important source of interstellar dust (Dwek 2004). 
The total mass of condensable elements (C, Mg, Si, S, Ca, 
Fe, Ni and associated O) ejected by a 30 M Q mass star is 
0.3-2 M leading to d/g ~ 0.01-0.07. Dust from red giants 



and AGB is less suitable to high velocity clouds. For com- 
parison, d/ g between <~ 2.5 10~ 3 and 6.6 10~ 3 were deter- 
mined for actively star-forming galaxies by Calzetti (2000) 
and for PG quasars by Haas et al (2003) adopting the Hilde- 
brand (1983) classical method. 

In Table 8 dust masses, M<j, calculated in the present 
work are compared with those calculated by the models of 
Farrah et al. (2002). Discrepancies within a factor < 3 are 
found for all the galaxies of the sample considering the er- 
rors, except for F10026+4949 and 18216+6418 where the 
discrepancies are greater than a factor of 10. The discrep- 
ancies depend on the fact that the dust masses were calcu- 
lated by Farrah et al. with the Hildebrand (1983) method. In 
18216+6418 we predict very low d/g, because large slabs of 
gas are deduced from the strong self-absorption of free-free 
radiation at long wavelengths (Fig. 7, last diagram). Adopt- 
ing a temperature T=1000 K and a density n=10 4 cm~ 3 
downstream of clouds in physical conditions corresponding 
to V s =500 km s" 1 and n =300 cm , the gas is optically 
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Table 8 Comparison of M d calculated in the present work with 
M d reported by Farrah et al. (2002) in HLIRGs 





Mi 

(10 8 M Q ) 


M d (Farrah et al.) 
(1O 8 M ) 


F00235+1024 


1.86 


1.95 


07380-2342 


17.8 


6.17 


F 10026+4949 


35. 


<5 


F 12509+3 122 


6.4 


1.58 


13279+3401 


0.22 




14026+4341 


2.0 


1.74 


F14218+3845 


1.4 


1.07 


F16124+3241 


6.1 


3.47 


EJ 1640+41 


1.9 


<1.58 


18216+6418 


0.3 


3.24 



1 adopting M gas = 10 11 M Q (Farrah et al. 2002) 

thick at v < 300 GHz for a layer > 242.2 pc thick (Oster- 
brock 1988), which is actually of the order of the calculated 
radius of the NLR (Table 6). 

6 Summary 

In this work we have modelled consistently the IR dust emis- 
sion and the bremsstrahlung. From the analysis of dust-to- 
gas ratios in galaxies with different IR luminosities, we con- 
clude that the high infrared luminosities generally depend 
on high amount of dust relative to gas close to the galaxy 
centre. This is revealed by the continuum emitted down- 
stream of relatively high velocity clouds. Relatively high 
d/g (0.04 by mass) were found in the LIRG clouds with 
V s = 300-500 kins" 1 , higher by a factor of - 100 than in 
the ISM (~ 4 10~ 4 ). 

Half of the LIRGs of the sample contain an AGN, which 
is found in almost all ULIRGs and in none of the HLIRGs. 

Interestingly, we could not demonstrate that higher Lir 
as those observed in HLIRGs correspond to higher dust- 
to-gas ratios. On the contrary, the highest d/g appear in 
LIRGs. Notice that the HLIRG types correspond to narrow- 
line galaxies, Seyfert 1, and QSO. It was found in previ- 
ous modelling (Contini et al. 2002) that d/g ratios in the 
narrow line region of these galaxies are lower by a factor 
of ~ 10 than for Seyfert 2 and LINERs. It seems that in 
some HLIRGs the high IR luminosities include the con- 
tribution of bremsstrahlung from cool gas. The relatively 
high bremsstrahlung and the large distances of the emitting 
clouds from the galaxy centre suggest that the HLIRGs of 
the Farrah et al sample were selected among large, massive 
objects. 
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